Operation of a Stark decelerator with optimum acceptance 
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With a Stark decelerator, beams of neutral polar molecules can be accelerated, guided at a con- 
stant velocity, or decelerated. The effectiveness of this process is determined by the 6D volume in 
phase space from which molecules are accepted by the Stark decelerator. Couplings between the 
longitudinal and transverse motion of the molecules in the decelerator can reduce this acceptance. 
These couplings are nearly absent when the decelerator operates such that only every third electric 
field stage is used for deceleration, while extra transverse focusing is provided by the intermediate 
stages. For many applications, the acceptance of a Stark decelerator in this so-called s = 3 mode 
significantly exceeds that of a decelerator in the conventionally used (s = 1) mode. This has been 
experimentally verified by passing a beam of OH radicals through a 2.6 meter long Stark decelerator. 
The experiments are in quantitative agreement with the results of trajectory calculations, and can 
qualitatively be explained with a simple model for the 6D acceptance. These results imply that the 
6D acceptance of a Stark decelerator in the s = 3 mode of operation approaches the optimum value, 
i.e. the value that is obtained when any couplings are neglected. 
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I. INTRODUCTION 

Since its introduction in 1999, the method of Stark 
deceleration has developed into an established method 
for taming molecular beams [l|, 0] ■ The Stark decelera- 
tion technique combines molecular beam technology with 
concepts from charged particle accelerator physics. In 
essence, a part of a beam of neutral molecules is selected 
and decelerated utilizing the force that polar molecules 
experience in inhomogeneous electric fields. A Stark de- 
celerator produces bunches of state-selected molecules 
with a computer-controlled velocity and with narrow ve- 
locity distributions. These beams are ideally suited for 
a variety of experiments in which the velocity of the 
molecules is an important parameter. Applications in- 
clude the use of slow molecular beams to enhance the 
interaction time in spectroscopic experiments [3, [j| and 
scattering studies as a function of the collision energy I5| . 
When the molecules are decelerated to a near standstill, 
they can be loaded and confined in traps [a, H]- This 
allows the observation of molecules in complete isolation 
for times up to several seconds, and enables the investi- 
gation of molecular properties in great detail [^, [3, |T3| ■ 

For many of these applications it is crucial that the 
number density of the decelerated packets of molecules 
is further increased. Higher densities of decelerated 
molecules will improve the statistics in metrology exper- 
iments and can be decisive, for instance, for the obser- 
vation of (in)elastic scattering or reactive collisions in 
crossed molecular beam experiments. Higher densities in 
the trap are also a prerequisite for the future applica- 
tion of cooling schemes like evaporative cooling, needed 
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to reach the regime of degenerate dipolar quantum gases 

The number density of decelerated molecules that can 
be reached at the exit of the decelerator is limited by 
the initial phase space density in the molecular beam 
and by the 6D phase space acceptance of the decelerator. 
The latter is defined as the volume in 6D phase space 
- the product of the volume in real space and in veloc- 
ity space - from which stable trajectories through the 
decelerator originate. In most Stark deceleration experi- 
ments to date, molecular beams with a low initial velocity 
are slowed down using decelerators with a rather limited 
number of electric field stages. Hence, molecular beams 
are typically released from a cooled pulsed valve using Xe 
or Kr as a carrier gas. The use of Xe or Kr and the cool- 
ing of the pulsed valve strongly enhances cluster forma- 
tion, however, and is generally regarded to be non-ideal 
for a molecular beam expansion. Moreover, these decel- 
erators are usually operated at large phase angles. The 
phase angle 0o determines the deceleration rate per elec- 
tric field stage, and ranges from 0° < (/)o < 90° for decel- 
eration, while acceleration occurs from —90° < 0o < 0°. 
While for increasing values of |(/)o| the deceleration rate 
gets larger, the longitudinal phase space acceptance gets 
smaller. 

An obvious route to higher number densities of decel- 
erated molecules is thus the use of seed gases of lower 
mass (preferably Ne or Ar) in room-temperature expan- 
sions and the use of low phase angles in the decelerator. 
Together this implies, however, that (much) longer Stark 
decelerators need to be constructed to compensate for, 
both, the higher initial velocity and the lower decelera- 
tion rate. It is not a priori clear whether one can actu- 
ally transport molecules through such long decelerators 
without significant losses. If only molecular trajectories 
along the molecular beam axis are considered, the length 
of the decelerator is inconsequential, as the deceleration 



process is subject to phase stability [l3,ll3|. In reality the 
molecules in the beam have off-axis position and veloc- 
ity components, and do not only experience forces in the 
longitudinal (forward) direction. The transverse electric 
field gradients in the decelerator drive the molecules back 
towards the molecular beam axis. The resulting trans- 
verse oscillatory motion is strongly coupled to the longi- 
tudinal motion and can result in a reduction of the 6D 
phase space acceptance of the decelerator [ij]. Numer- 
ical simulations indicate that this coupled motion does 
not affect the overall performance of the relatively short 
Stark decelerators that have been used and operated at 
high phase angles thus far, but that it can severely af- 
fect the performance of longer decelerators that are op- 
erated at low phase angles. In the ideal case, the lon- 
gitudinal and transverse motions in the Stark decelera- 
tor are completely uncoupled. This can be achieved by 
constructing decelerators with dedicated, spatially sepa- 
rated, elements for focusing and deceleration ^, il6iLas 
is common practice in charged particle accelerators [17| . 
The required electrode geometries make the decelerator 
rather complex, however [l^ [l3| . 

In this paper we exploit that in a Stark decelerator 
with the original electrode geometry [l[ , the coupling be- 
tween the longitudinal and transverse motion can be sig- 
nificantly reduced when the decelerator is operated in 
the so-called s = 3 mode [1^. In this mode, only every 
third electric field stage is used for deceleration, while 
extra transverse focusing is provided by the intermediate 
stages. We demonstrate and quantify that for many ap- 
plications, the acceptance of the Stark decelerator in the 
s = 3 mode significantly exceeds that of a decelerator in 
the conventionally used (s = 1) mode of operation. The 
improved performance of the s = 3 operation mode was 
demonstrated earlier for guiding at a constant velocity 
{(j)a = 0°) [14] and for deceleration in a relatively short 
decelerator at high phase angles [16|, but was not quan- 
tified thus far. Here we present experiments in which 
a beam of OH radicals passes through a 2.6 meter long 
Stark decelerator consisting of 316 electric field stages. In 
this machine, the OH radicals can be detected either af- 
ter 103 or after 316 electric field stages, enabhng a direct 
comparison between the s = 1 and s — 3 modes of oper- 
ation under otherwise identical conditions, in particular 
using the same phase angle 0o- The experimental results 
are in quantitative agreement with the results of trajec- 
tory calculations, and can qualitatively be explained with 
a simple model for the 6D acceptance. 



II. EXPERIMENT 

A. Experimental set-up 

The experimental setup is schematically shown in Fig- 
ure [H A pulsed beam of OH radicals is produced via 
ArF-laser (193 nm) dissociation of HNO3 seeded in an 
inert gas. The dissociation is carried out inside a quartz 



capillary that is mounted on the orifice of a pulsed valve 
(General Valve, Series 99). The experiment runs at a rep- 
etition frequency of 10 Hz. During the supersonic expan- 
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FIG. 1: Scheme of the experimental setup. A pulsed beam 
of OH radicals is produced via photodissociation of HNO3 
seeded in Xe, Kr, or Ar. The beam of OH radicals passes 
through a 2.6 meter long Stark decelerator that consists of 
three modules of ~ 100 stages each. The OH radicals can be 
state-selectively detected using a laser induced fluorescence 
scheme at the end of the decelerator, and in the region be- 
tween the first two modules. In the top inset, a photograph 
of a decelerator module is shown. 

sion, the majority of the OH radicals cool to the lowest 
rotational and vibrational level of the X ^^zl2 spin-orbit 
manifold of the electronic ground state. This population 
is equally distributed over the two A-doublet components 
of the J — 3/2 level. Only OH molecules in the upper A- 
doublet component are low-field-seeking, and are of rele- 
vance to the experiments discussed here. This component 
splits into a Af/il — —3/4 and a Afjfi = —9/4 compo- 
nent in an electric field. Molecules in the Mj^ = —9/4 
component experience a Stark shift that is a factor of 
three larger than the Stark shift experienced by molecules 
in the AfjJ7 = —3/4 component. 

After passage through a skimmer with a 2 mm diame- 
ter opening, the molecular beam enters the differentially 
pumped decelerator chamber. The skimmer is mounted 
on a compact gate valve [20| , enabling the venting of the 
source chamber while keeping the decelerator chamber 
under vacuum. The beam enters the Stark decelerator 
60 mm from the nozzle orifice. The Stark decelerator 
consists of three modules that are mechanically and elec- 
trically decoupled from each other. The first two mod- 
ules consist of 104 electrode pairs (i.e. 103 electric field 
stages) each, whereas the last module contains 109 elec- 
trode pairs. These electrode pairs consist of two parallel 
4.5 mm diameter stainless steel electrodes that are placed 
symmetrically around the molecular beam axis, provid- 
ing a 3 mm gap for the molecular beam to pass through. 
Adjacent electrode pairs are alternately horizontally and 
vertically oriented, such that a 3 x 3 mni^ opening area 



remains for the molecular beam. The electrodes of all 
horizontal (vertical) pairs within each module are elec- 
trically connected and switched simultaneously to high 
voltage. The center-to-center distance (L) of electrodes 
of adjacent pairs is 8.25 mm, and the three modules are 
carefully aligned to also maintain this distance between 
the electrode pairs of adjacent modules. The first and 
last eight electrode pairs of each module are mounted 
on conically shaped rods, as shown enlarged in Figure 
[TJ This design provides excellent optical access for fluo- 
rescence collection in between adjacent decelerator mod- 
ules. It also allows the exits of two Stark decelerators in a 
crossed beam arrangement to be brought close together. 

The electric field in the decelerator is switched back 
and forth between two different configurations, that are 
schematically shown in Figure [21 In each configuration, 
the opposing electrodes of every other electrode pair are 
at ±20 kV, while the remaining electrodes are grounded. 
Switching between the two configurations is performed 
using fast air-cooled high voltage switches (Behlke HTS 
30I-03-GSM). To minimize the power dissipation per 
switch, each module is connected to four separate high 
voltage switches. Each switch is connected to its power 
supply via a 0.5 fjF capacitor bank, limiting the voltage 
drop during a deceleration cycle to less than 5 %. For 
the conventional mode of operation of a Stark decelera- 
tor, the s = 1 mode, the voltages are switched every time 
the molecules approach the pair of electrodes that are on 
high voltage [l3|- In this case, the molecules are simul- 
taneously decelerated and transversally focussed in every 
electric field stage. When the decelerator is operated in 
the s = 3 mode, the voltages are switched only after every 
third stage. In this case, only every third stage is used 
for combined deceleration and transverse focusing, while 
the intermediate stages provide additional focusing. 

The OH radicals can be state-selectively detected using 
an off-resonant Laser Induced Fluorescence (LIF) detec- 
tion scheme at two different positions along the beam 
line. The first detection zone is located between the first 
two modules and the second one is 18 mm downstream 
from the last module, enabling the detection of OH rad- 
icals after 103 or 316 electric field stages, respectively. 
The 282 nm radiation of a pulsed dye laser crosses the 
molecular beam in either one of the detection regions at 
right angles, and saturates the (spectroscopically not re- 
solved) Q2i(l) and Qi(l) transitions of the A^T,^,v — 
1 <— X ^11^/2, V = band. The fluorescence occurs pre- 
dominantly on the yl^I]+,w = 1 ^ X^H, u = 1 transi- 
tion around 313 nm. Stray light from the laser is min- 
imized by passing the laser beam through light baffles 
between the entrance and exit windows, and by optical 
filtering in front of the photomultiplier tube (PMT) . 

In the experiments the seed gases Xe, Kr, or Ar are 
used. The mean initial beam velocities and velocity dis- 
tributions are deduced from time-of-flight (TOF) mea- 
surements of the OH beam when all electrode pairs of 
the decelerator are simultaneously charged to ± 7 kV. 
The mean velocities for the beam are 350 m/s, 430 m/s. 
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FIG. 2: Schematic representation of the two electric field con- 
figurations that are used in the deceleration process, together 
with the potential energy W for an OH molecule along the 
molecular beam axis. By switching between the two field con- 
figurations when the molecules are at the positions indicated 
by the vertical dashed lines, an amount of kinetic energy AW 
is removed from the molecules. In the conventional (s — 1) 
mode of operation, each electric field stage is used simultane- 
ously for deceleration (D) and focusing in alternating trans- 
verse directions {Fx,Fy). In the s = 3 mode of operation, 
only every third stage is used for combined deceleration and 
transverse focusing, while the intermediate stages provide ad- 
ditional focusing. 



and 590 m/s for Xe, Kr, and Ar, respectively, with a 
full width at half maximum (FWHM) velocity spread of 
about 15-20 % for all seed gases. 

In Figure [3] the intensity of the LIF signal of a beam 
of OH (J — 3/2) radicals seeded in Xe is shown as a 
function of time after firing the dissociation laser using 
different deceleration sequences. The OH radicals are de- 
tected using the second LIF detection unit, and the Stark 
decelerator is operated in the s = 3 mode. In curve (b) 
the TOF profile is shown that is obtained when the de- 
celerator is operated at a phase angle of (po = 0°, corre- 
sponding to guiding a packet of OH radicals at a constant 
velocity. A packet of OH radicals with a mean velocity 
of 350 m/s is selected, transported through the 2.6 me- 
ter long decelerator, and arrives in the detection region 
some 7.6 ms after its production, with a FWHM of the 
arrival time distribution of 25 /iS. The measurements 
shown in curve (a) are obtained with the decelerator op- 
erating at a phase angle of 0o = —50°, accelerating a 
packet of OH radicals from an initial velocity of 350 m/s 
to a final velocity of 544 m/s. There is no signature of 
the part of the molecular beam that is not accelerated. 
This is also expected as the electrodes of the decelera- 
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FIG. 3: Time-of-flight profiles of OH radicals, recorded at the 
exit of the 316 stage Stark decelerator. The Stark decelerator 
is operated at the s — 3 mode, and accelerates (curve (a)), 
guides (curve (6)), or decelerates (curve (c)) a packet of OH 
radicals with an initial velocity of 350 m/s. The TOF pro- 
files that result from simulations of the experiment are shown 
underneath the experimental profiles. 



tor are switched to ground when the accelerated packet 
exits the decelerator, about 1.5 ms before the remainder 
of the beam pulse would arrive in the detection region. 
Curve (c) shows the TOF profile that is obtained when 
the decelerator is operated at a phase angle of cpo = 29° 
to decelerate a packet of OH radicals from 350 m/s to 
150 m/s. The decelerated molecules exit the decelera- 
tor about 10.7 ms after production, about 3 ms after the 
arrival of the undecelerated part of the beam. 

The experimental TOF profiles are in excellent agree- 
ment with the profiles that result from three dimensional 
trajectory simulations of the experiment that are shown 
underneath the experimental profiles. In these and in all 
subsequent simulations, the individual contributions of 
the Mjft — —3/4 and the Mjft = —9/4 components to 
the LIF signal intensity are taken into account. 



B. Comparing the s = 1 and s — 3 modes of 
operation 

In Figure S] two series of TOF profiles are shown that 
allow a direct comparison between the performance of 
a Stark decelerator in the s = 1 and s = 3 operation 
mode under otherwise identical conditions. In both se- 
ries, the Stark decelerator is programmed to accelerate, 
guide or decelerate a packet of OH radicals with a mean 
initial velocity of 350 m/s to a final velocity that is in 
the 70 - 600 m/s range (-65° < (po < 32°). Only that 
part of each TOF profile that contains the signature of 
the packet at the final velocity is shown. In the left and 
right panels the series of profiles are shown that are ob- 
tained when the decelerator is operated using the s = 1 
and the s = 3 mode, detecting the OH radicals after 



103 and 316 stages, respectively. The (almost) factor of 
three difference in the number of stages results in (al- 
most) identical phase angles for the Stark decelerator to 
produce a given final velocity in both series. The phase 
angle that is used, and the final velocity of the packet, is 
indicated for selected profiles in both panels. To enable a 
direct comparison between the two modes of operation, 
both series are plotted on the same vertical scale. For 
this, the relative detection efficiency in the two LIF zones 
is experimentally determined by measuring OH radicals 
at both detection locations when the Stark decelerator 
is operated at s — 3,(f>o — 0°. The overall scaling factor 
that is thus determined is correct if we assume that for 
the s = 3,(j)o = 0° mode of operation the number density 
of the packet of OH radicals does not decline when pro- 
gressing from the first to the second LIF zone, and that 
the relative detection efficiency is independent of the ve- 
locity of the OH radicals. The validity of both assump- 
tions is checked experimentally and verified by numerical 
simulations. 

When the decelerator is operated using s = 1, the sig- 
nal intensity for 00 = 0° is about a factor of two-and-a- 
half lower than the signal intensity for s ~ 3,0o = 0°. 
The signal intensity for the s = 1 appears rather constant 
for the different values of (pQ that are used. For the s — 3 
mode of operation it is observed that the signal intensity 
for (f>Q y^ 0° gradually reduces, following the reduced ac- 
ceptance of the decelerator for increasing absolute values 
of the phase angle. For final velocities below 150 m/s, in- 
dicated by the dashed line in Figure [4j a sharp reduction 
of the signal intensity is observed even though the change 
in phase angle is only very small. This reduction is due to 
excessive transverse focusing for low velocities, and will 
be discussed in more detail in section [HDJ The signal in- 
tensities for the s = 1 and s = 3 modes of operation are 
about equal for 0o — —65°. Both series of TOF profiles 
are in excellent agreement with profiles that result from 
simulations of the experiment that are shown underneath 
the experimental profiles. It is noted that the low field 
seeking Mjft — —3/4 component only contributes to the 
TOF profiles for -30° < (/)o < 30°. 

The gain of the s — 3 mode with respect to the s = 1 
mode of operation, defined as the ratio of the maximum 
signal intensities of the decelerated packets at a given 
final velocity, is shown as a function of the final velocity 
in Figure [5] For this, the data presented in Figure 2] are 
used. For selected velocities, the phase angle that is used 
for s — 3 operation is indicated; the phase angle that 
is used for s = 1 operation differs only slightly from this 
value. It is seen that for (po = 0°, i.e. at a final velocity of 
350 m/s, the gain is 2.4, consistent with previous studies 
(la , [l9| . When the beam is accelerated to 450 m/s (0o — 
—21.7°), a gain up to 3 is observed. For higher velocities 
the gain gradually reduces and reaches 1.2 for 600 m/s 
((/)o — —63.2°). A gain up to 4.2 is observed when the 
beam is decelerated to 200 m/s (00 — 22.3°). Below 
150 m/s (00 — 27.0°), the gain drops fast and reaches 
1.0 for a final velocity of 120 m/s. 
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FIG. 4: Time-of-flight profiles of OH radicals that exit the Stark decelerator using the s = 1 and s — 3 mode of operation. 
The OH radicals are detected after 103 and after 316 electric field stages for s = 1 and s — 3, respectively. The measurements 
are recorded under otherwise identical conditions, and are shown on the same vertical scale. The beam of OH radicals has a 
mean initial velocity of 350 m/s. The mean final velocity of the molecules and the phase angle used, are indicated for selected 
profiles. 



C. The s = 1 mode of operation at low phase angles 
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FIG. 5: The ratio of the maximum signal intensities (squares, 
connected by straight line segments) for s — 3 operation ver- 
sus s — 1 operation as a function of the final velocity. The two 
additional data points (stars) apply to a bimodal operation 
of the decelerator (see text for details). 



In this section we want to address the question 
whether, for a fixed initial and final velocity, the number 
density of decelerated molecules can also be increased by 
using lower phase angles in the s = 1 mode of operation. 
For this, we compare the deceleration with 103 stages at a 
certain phase-angle to the deceleration with 316 stages at 
about one-third of this phase-angle. The data for the de- 
celeration at the s = 1 mode with 103 deceleration stages, 
have already been shown and discussed in the previous 
section. The complimentary data for the deceleration at 
the s — 1 mode with 316 deceleration stages have been 
measured as well; a beam of OH radicals with the same 
initial velocity of 350 m/s is decelerated or accelerated 
to the same final velocities between 70 m/s and 600 m/s. 
The phase-angles used for these measurements are very 
low, ranging from —21.7° to 10.9°. Again, the signal in- 
tensity that is observed with guiding at the s — 3 mode, 
is used to calibrate the measurements with 316 stages 
relative to the ones with 103 stages. 

For the operation in the s — 1 mode, the gain in us- 
ing 316 stages compared to using 103 stages is shown 
as a function of the final velocity in the lower curve in 
Figure [HI This gain is determined as the ratio of the 
signal intensities of the decelerated packets at a given 
final velocity, and lies between 0.5 and 0.7 throughout. 
Since this gain stays smaller than one in the considered 
velocity range, it is evident that for the chosen parame- 



ters the number density of decelerated molecules can not 
be increased by the use of lower phase-angles. This per- 
haps counterintuitive finding is explained by the presence 
of inherent instabilities in the s = 1 mode, which more 
strongly manifest themselves during the increased time 
spent in the decelerator. 
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FIG. 6; Lower curve: The ratio of the signal intensity for 
s = 1 operation of a decelerator with 316 versus 103 stages 
as a function of the final velocity. Upper curve: The ratio of 
the signal intensity using a 316 stage decelerator operating on 
the s = 3 versus the s = 1 mode. 

We can now also address the question whether, given 
a Stark decelerator of a certain length, the maximum 
number density of decelerated molecules is obtained using 
the s = 3 mode at a certain phase-angle or using the s = 1 
mode at about one-third of this phase-angle. For the 316 
stage decelerator, the resulting gain of the s — 3 versus 
the s = 1 mode directly follows from the curve shown in 
Figure [5] and the lower curve in Figure [SI and is shown as 
the upper curve in the Figure [6] Operation on the s — 3 
mode is seen to be about a factor five better than on the 
s = 1 mode, provided the final velocity is larger than the 
threshold velocity. 



D. Excessive focusing at low^ velocities 

The rather abrupt decrease in the number density of 
decelerated molecules for velocities below 150 m/s as 
shown in the right hand panel of Figure 2] can quali- 
tatively be understood as follows. During their flight 
through the decelerator, molecules are alternatingly fo- 
cused in each transverse direction. When the focus- 
ing force acts in one transverse direction, molecules ex- 
perience to a good approximation no focusing or de- 
focusing force in the orthogonal transverse direction (the 
molecules actually experience a small defocusing force 
in the orthogonal direction). As long as the charac- 
teristic wavelength A of the transverse oscillatory mo- 



tion is much larger than the periodicity of the transverse 
focusing force, molecules will follow stable trajectories 
through the decelerator. The wavelength A is given by 
A = (vz) 27r/r2j,, where (vz) and fly are the mean longi- 
tudinal velocity and the mean transverse oscillation fre- 
quency of the packet of molecules, respectively. The pe- 
riodicity of the focusing force is given by 2sL, where L is 
the center-to-center distance of adjacent electrode pairs. 
For high velocities, therefore, stable trajectories are ex- 
pected. For low velocities, however, A becomes ever closer 
to 2sL, and molecules will get more tightly transversely 
focused. The molecular trajectories will then exhibit ever 
larger deviations from the molecular beam axis, and the 
molecules will eventually crash onto the electrodes. For 
a given electric field distribution in a Stark decelerator, 
the resulting loss of molecules is thus expected to strongly 
depend on the longitudinal velocity of the molecules. 

In Figure [71 the maximum signal intensity of decel- 
erated packets of OH radicals is shown at the exit of 
the 316 stage decelerator, operating in the s = 3 mode. 
Beams of OH radicals with three different initial veloci- 
ties have been used, and the signal is shown as a function 
of the final velocity (upper panel) or as a function of the 
phase-angle (lower panel). When Xe or Kr are used as 
carrier gas, the Stark decelerator is programmed to se- 
lect a packet of molecules with an initial velocity that 
is identical to the mean velocity of the molecular beam, 
i.e. 350 m/s or 430 m/s, respectively. For Ar, a ve- 
locity of 520 m/s is selected from the slow tail of the 
velocity distribution of the beam. The phase angle (j)o is 
varied to decelerate the selected packet of molecules to 
final velocities down to 100 m/s. The series of measure- 
ments for each seed gas arc normalized to the data point 
that corresponds to 0o = 0°. The thresholds are found 
at a velocity of about 150 m/s in each series, and are 
indicated by the vertical dashed lines. This velocity is 
reached when 0o = 27°, (j)o — 43°, and 4>o — 67° for Xe, 
Kr, and Ar seeded beams, respectively. The value of the 
threshold velocity appears independent from the phase 
angle (J)q of the decelerator, consistent with the qualita- 
tive picture described above. The intensities that result 
from numerical trajectory simulations of the experiments 
are shown as solid lines in Figure [71 Excellent agreement 
is obtained with the experiments and in particular the 
threshold behavior of the signal intensity is reproduced 
well. 

The threshold velocity below which losses due to ex- 
cessive transverse focusing occurs, can be approximately 
determined as described in Appendix B. The model de- 
scribed there predicts that no stable trajectories exist in 
the s ~ 3 mode of operation when A < (2.6 • 2sL). For 
the present Stark decelerator, operating on s = 3, the 
periodicity of the transverse focusing force 2sL is equal 
to 49.5 mm. The mean transverse oscillation frequency 
fiy follows from the time-averaged transverse force (see 
Appendix A). For s = 3, fiy is rather independent from 
the phase angle (j)o, and a threshold velocity of about 
92 m/s is found, consistent with the experimental find- 
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FIG. 7: Maximum signal intensity of decelerated packets of 
OH radicals as a function of the final velocity (upper panel) 
and as a function of (j>o (lower panel), using a 316 stage de- 
celerator operating at the s = 3 mode. Beams of OH radi- 
cals with three different initial velocities, produced by seeding 
in Xe, Kr, or Ar, have been used. The intensities that are 
obtained from numerical simulations of the experiment are 
shown as solid lines. 



ings. When the decelerator is operated in the s — 1 
mode, a threshold velocity is not so clearly defined, but 
similar losses due to excessive transverse focusing occur 
for velocities below 30 m/s. 

The rather high threshold velocity for s — 3 does not 
severely affect experiments in which Stark-decelerated 
beams are used for high resolution spectroscopy and col- 
lision studies, or in which the decelerated beams are in- 
jected into molecular storage rings or synchrotrons. It 
does affect, however, experiments in which lower final 
velocities are required, e.g. trap loading experiments. 
There are several approaches to yet produce decelerated 
packets at a velocity below the threshold velocity with 
decelerators that are intended to operate at s = 3. An 
electric field geometry for the last section of the deceler- 
ator can be designed that permits a gradual reduction of 
the transverse focusing strength. This can be achieved 
by a dedicated electrode geometry and/or by a sequential 



reduction of the voltage that is applied to the electrodes. 
It is noted that similar strategies have already been im- 
plemented in tra ppin g experiments using decelerators in 
the s = 1 mode [2H. An alternative approach is to de- 
velop an electrode geometry for the last segment of the 
decelerator that allows the confinement of molecules in 
a genuine traveling potential well. When the velocity of 
this well is gradually reduced, the packet of molecules can 
be transferred from the threshold velocity to lower veloc- 
ities without loss. The trapping of molecules in genuine 
traveling potential wells has already been demonstrated 
using optical fields [231 and using electric fields above a 
micro-structured electrode array fc23j |. 

Within the possibilities of the present experimental ar- 
rangement, decelerated packets with low final velocities 
can be produced by changing over from the s = 3 mode 
to the s = 1 mode before the threshold velocity has been 
reached. The number of molecules that exit the decel- 
erator at velocities below the threshold velocity strongly 
depends on the details of the change-over, i.e., the veloc- 
ity after which s = 1 operation is used, and the phase 
angles that are used before and after the transition. The 
influence of the choice of these parameters on the num- 
ber of molecules that exit the decelerator has experimen- 
tally been studied, decelerating a packet of OH radicals 
with an initial velocity of 350 m/s. The velocity and 
position of the packet in the decelerator at which the 
change-over from s = 3 to s = 1 is made is systemati- 
cally varied for the target (flnal) velocities of 100 m/s and 
50 m/s. For both velocities the maximum signal intensity 
is observed when the transition to s = 1 is made when 
the first ~ 300 stages are operated at s = 3, 0o ~ 26° 
and when the molecular packet has reached a velocity of 
170 m/s. The remaining 12 and 15 stages are then used 
a,t s = l,(j>o = 43.2° and s = l,0o = 48.3° to produce 
the final velocities of 100 m/s and 50 m/s, respectively. 
The gain of this bimodal operation of the Stark deceler- 
ator with respect to s = 1 operation is shown in Figure 
[5] as separate data points. It is observed that for a fi- 
nal velocity of 100 m/s the gain is about a factor of 2, 
and is close to one for a final velocity of 50 m/s. These 
measurements demonstrate that low final velocities can 
be produced with Stark decelerators that are designed 
to operate at s = 3, although the efficiency approaches 
the efficiency of conventional Stark decelerators when the 
target velocity is much lower than the s = 3 threshold 
velocity. 



III. NUMERICAL TRAJECTORY 
CALCULATIONS 

The measurements presented thus far demonstrate the 
performance of the s = 1 and s = 3 modes of operation 
of Stark decelerators in the (limited) range of parame- 
ters that is accessible to the experiment. In this section, 
both modes of operation are studied in a wider param- 
eter range using numerical trajectory simulations. The 



electrode geometry that is used in the simulations is the 
same as used in the experiments, but the decelerator is 
allowed to have an arbitrary length. 

Trajectories of OH radicals through the decelerator are 
numerically calculated as a function of the phase angle 
00 for both the s = 1 and the s = 3 mode of opera- 
tion. In these simulations, a large number of molecules is 
homogeneously distributed at the entrance of the Stark 
decelerator over a block in 6D phase space. This block 
has a dimension of 20 mm x 90 m/s in the longitudinal 
direction, and a dimension of 4 mm x 25 m/s in each 
transverse direction. The molecular distribution has a 
mean forward velocity of 550 m/s at the entrance of the 
decelerator, and is decelerated or accelerated to a final ve- 
locity of 180 m/s or 755 m/s, corresponding to a change 
in the kinetic energy of 90 %. These values are arbi- 
trary and can be chosen without loss of generality, as the 
phase-space acceptance of a Stark decelerator is in princi- 
ple independent of the absolute initial and final velocity. 
The rather high final velocity of 180 m/s for decelera- 
tion is chosen to stay away from the velocities for which 
excessive transverse focusing occurs, as discussed in sec- 
tion llTDl Decelerators containing 781, 388, 256, 190, 151, 
126, 111, 102, and 99 electric field stages are simulated 
that are operated using \(j)o\ = 10°, 20°, 30°, 40°, 50°, 
60°, 70°, 80°, and 90°, respectively. These numbers ap- 
ply to s = 1; for simulations that apply to s = 3 the num- 
ber of stages is three times as large. For both values of s, 
additional simulations were performed for 0o = 0° using 
a 2500 stages long decelerator. In each simulation, a suf- 
ficient number (5.000.000 for s = 1 and 500.000 for s = 3) 
of molecules is generated to obtain good statistics. The 
number of molecules that are within the phase-space dis- 
tributions of the decelerated packet are counted, and the 
corresponding 6D volume in phase-space is calculated. 
In the lower part of Figure [51 the simulated longitudinal 
phase-space distributions for (j)o = —10° are shown both 
for s = 1 and for s = 3, together with the separatrices 
that follow from the ID model for phase stability [19J . 
These distributions are representative for the distribu- 
tions at low phase angles in general, and are shown here 
to exemplify the simulation method only. The phase- 
space distribution for 00 = — 10°,s = 1 is highly struc- 
tured with alternating stable and unstable regions. In 
the distribution for 00 — — 10°,s = 3 no clear structure 
is present. The area within the longitudinal separatrix ~ 
the longitudinal acceptance - is a factor v3 smaller for 
s = 3 than for s ~ 1 [ig |. 

In the upper panel of Figure [8l the resulting simulated 
6D phase-space acceptance is shown as a function of 0o 
for both s = 1 and s — 3. For s = 1 the phase-space 
acceptance is maximum for 0o = 0°, and drops signifi- 
cantly for 00 7^ 0°. It has a minimum at |0o| = 10°, is 
rather constant for 20° < |0o| < 50°, and drops again 
for 1 00 1 > 50°, consistent with the experiments in sec- 
tion [ircl Note that the Stark decelerator has a nonzero 
acceptance for |0o| = 90°, as discussed before [ij. It 
is interesting to note that the phase-space acceptance is 
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FIG. 8: 6D phase space acceptance of a decelerator as a func- 
tion of tlie phase- angle 0o, resulting from numerical trajectory 
calculations (squares connected with solid lines for s = 3; tri- 
angles connected with dashed lines for s = 1), together with 
the model predictions (solid line for s = 3; dashed line for 
s = 1). In the inset, the s = 1 data is shown enlarged. In the 
lower part, the longitudinal phase-space distributions that re- 
sult from the simulations are shown for 0o = —10°, both for 
s = \ and s = 3. 



not symmetric around 0o = 0°. The phase-space accep- 
tance for 00 < is larger than the acceptance for the 
same positive phase-angle; molecules spend less time in 
the decelerator when they are accelerated then when they 
are decelerated, reducing the loss due to instabilities. 

The calculated acceptance for s — i shows a rather 
different trend, and is generally larger for smaller values 
of |0o|. Although the longitudinal phase-space accep- 
tance for s = 3 is a factor -s/S smaller than for s = 1, 
the 6D phase-space acceptance exceeds the acceptance 
for s = 1 for 0° < |0o| < 70°. When |0o| = 70°, the 
acceptance in both modes of operation is equal, consis- 
tent with the experimental findings discussed in section 
IIIBI For |0o| > 70° the acceptance for s = 1 is slightly 
larger than the acceptance for s = 3. Although less pro- 
nounced, the 1 00 1 dependence of the acceptance for s = 3 
is again asymmetric around 0o = 0°. 

It is interesting to compare the calculated phase-space 
acceptance with the phase-space acceptance that is ex- 
pected for a Stark decelerator in which instabilities are 
absent. From this comparison one can quantify the pres- 



ence and severeness of instabilities in the s = 1 and s = 3 
mode of operation. The longitudinal phase-space accep- 
tance of a decelerator is easily calculated and is, as men- 
tioned above, simply given by the area within the sepa- 
ratrix [ij, [l^ . The equations of motion that govern the 
transverse trajectories of molecules through the deceler- 
ator contain time-dependent forces, for which in general 
no simple analytical solutions exist. These equations, 
together with the equation for the longitudinal motion, 
can be used to estimate the volume in phase space from 
which stable trajectories can originate. This procedure 
is outlined in detail in Appendix A, and the resulting 
6D phase-space acceptance is shown as a function of (J)q 
for both s — I and s = 3 in Figure [H It is seen that for 
s = 1 the phase-space acceptance predicted by the model 
deviates significantly from the calculated acceptance; for 
I'/'ol < 20° the deviation is at least an order of magni- 
tude. For larger angles the discrepancy gets less, and 
both curves cross around |0o| = 70°. From this compar- 
ison it is once more evident that instabilities are present 
when the decelerator operates in the s = 1 mode [ij], 
and that these instabilities severely limit the obtainable 
acceptance. 

For the s = 3 mode, the acceptance predicted by the 
model reproduces the calculated acceptance much better. 
The agreement for |(/)o| > 40° is good, and the deviations 
are in the 5-20% range for 20° < |0o| < 40° and about 30 
% for |0o| < 20°. These minor deviations can be taken as 
an indication for the presence of small instable regions, 
as have indeed been observed close to the separatrix for 
00 = 0° [lj|. The overall agreement, however, demon- 
strates that the 6D acceptance of a Stark decelerator in 
the s = 3 mode of operation approaches the optimum 
value, i.e. the value that is predicted from the model 
that neglects any instabilities. 



IV. CONCLUSIONS 

The studies presented in this paper address the ques- 
tion how one can get the highest number density of decel- 
erated molecules with a certain velocity at the exit of a 
decelerator. Rather than discussing a variety of electrode 
geometries that one might use to decelerate a beam of po- 
lar molecules, these studies focus on a Stark decelerator 
in the conventional, experimentally proven design. This 
decelerator can run at different phase-angles and oper- 
ate in different modes, and can be built with a variable 
length. The number density of accelerated and decel- 
erated OH radicals has been experimentally studied as a 
function of these three parameters. Quantitative compar- 
isons of these number densities, obtained using Stark de- 
celerators with different parameter sets, have been made. 
The measurements have been substantiated by numeri- 
cal simulations, from which comparisons for a much wider 
range of parameters can be made. These studies provide 
quantitative arguments for the design criteria of Stark 
decelerators for specific applications. 



Based on the one-dimensional description of a Stark 
decelerator, one would expect more molecules at the end 
of the decelerator for longer decelerators that run at lower 
phase-angles. This description neglects the coupling be- 
tween the longitudinal and transverse motion, however, 
which limits the actual 6D acceptance of a decelerator. A 
first important conclusion from the present study is that, 
for a decelerator operating in the s = 1 mode, a strat- 
egy to optimize the number of decelerated molecules by 
using low phase angles and a large number of decelera- 
tion stages is only of limited use. There is a maximum 
of the 6D acceptance for a phase- angle of around 50°, 
and the optimum number of molecules is obtained when 
the length of the decelerator (for a given initial and final 
velocity) is adjusted such that this phase-angle can be 
used. A decelerator of 150 stages that is operated at 50°, 
for instance, produces more decelerated molecules than 
a decelerator of 250 stages that runs at 30° . 

A second important conclusion is that a decelerator 
that operates in the s = 3 mode outperforms a decelera- 
tor in the s = 1 mode in almost all cases. In the s = 3 
mode, coupling between the longitudinal and transverse 
motion is nearly absent, and lower phase angles always 
result in a larger acceptance. For small phase angles, a 
gain up to a factor of ten can be obtained. The gain 
depends strongly on the phase angle that is used and for 
phase angles above 70°, the acceptance for the two modes 
of operation is very similar. An intrinsic disadvantage of 
the s — 3 mode is that there are large losses for final ve- 
locities below around 150 m/s. Lower velocities can still 
be produced, however, and different schemes have been 
discussed and demonstrated for this. 



A third important conclusion is that the acceptance 
of a Stark decelerator operating in the s = 3 mode ap- 
proaches the optimum value, i.e. a decelerator operating 
in this mode works as good as it can be expected to get. 
This conclusion is based on a comparison between the 
outcome of numerical trajectory calculations and the 6D 
acceptance that is derived from a model. In this model, 
couplings between the longitudinal and transverse mo- 
tion are neglected. 

To make use of the advantages that the s — 3 mode of 
operation offers, a considerably longer Stark decelerator 
is needed than for the s = 1 mode. This indeed requires 
more electrode pairs and a longer vacuum chamber, but 
it should be realized that there is no additional require- 
ment on the high voltage electronics. Compared to the 
decelerators that have been commonly used so far (s = 1, 
« 100 stages, 4>q — 50 — 60°) a five times longer version 
operating in the s — 3 mode at somewhat lower phase 
angles will typically result in a factor five higher num- 
ber density at the exit. Moreover, this gain in number 
density is accompanied by a reduction in the longitudinal 
translational temperature. 
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Appendix A 

In previous work we have studied the transverse stabil- 
ity in a Stark decelerator [Ij]. Here we adapt the model 
for the transverse motion of molecules through a Stark 
decelerator, with the goal to derive the 6D phase-space 
acceptance as a function of the phase angle (/jq that can be 
expected when the presence of instabilities is neglected. 

In the description of the motion of the OH radicals 
through the decelerator, the z coordinate describes the 
position of the molecule along the molecular beam axis, 
while X and y are the transverse coordinates. The forces 
in the x and the y direction are assumed to be uncoupled 
from each other and identical. The alternating focussing 
in either one of the transverse directions (say y) is repre- 
sented by a continuously acting transverse force Fy (0, y) 
that depends on the phase 4> of the molecule and on y 

t+2T 

Fy{<P,y) = ^ ( F.yiy{t'),z{t'),t')dt' 



2T 

1 
2lL 



(0+2sTr)L/7i 



Fy{y,z)dz, (1) 



0L/71 



where 2T is the time during which the synchronous 
molecule travels a distance 2sL. 

To a good approximation, the transverse force Fy is lin- 
ear in the displacement y from the molecular beam axis. 
The strength of the transverse force can be expressed 
in terms of a frequency a;y(0)/27r, referred to hereafter 
as the natural transverse oscillation frequency, using the 
relation: 



Fv{(t>,y) ^-rnujl{(t>)y, 



(2) 



where m is the mass of the OH radical. In Figure [9l the 
natural transverse oscillation frequency is shown for an 
OH (A 2H3/2, J = 3/2, Mjfl = -9/4) radical as a func- 
tion of its phase </>. For s = 1 , the natural transverse os- 
cillation frequency has a strong dependence on the phase 
(j). For molecules close to (/i = 0°, the transverse fre- 
quency is very low and focussing forces are almost absent. 
For s = 3, the natural transverse oscillation frequency is 
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FIG. 9: Natural transverse oscillation frequency L0y/2n for an 
OH (X^Hg/a, J = 3/2, Mjfi = -9/4) radical as a function 
of its phase <j> in a Stark decelerator, for the operation modes 
s = 1 and s = 3. 



rather independent from the phase (j). The intermediate 
electric field stage focuses all molecules equally efficient, 
effectively decoupling the transverse focusing properties 
of the decelerator from the longitudinal phase 4>. 

The transverse phase-space acceptance for a given 
mode of operation of the decelerator is evaluated as 
schematically illustrated in Figure [101 The decelerator 
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FIG. 10: Schematic representation of the method used to cal- 
culate the transverse phase-space acceptance, (a) The tra- 
jectory of a molecule in longitudinal phase-space, shown as a 
dashed line, (b) The time dependence of the natural trans- 
verse oscillation frequency is constructed from the time de- 
pendence of the phase <^(t), and the phase dependence of the 
natural transverse oscillation frequency uJy{(f>)/2'K. (c) The 
time-averaged transverse focusing force results in elliptical or- 
bits in transverse phase-space. 

is operated at (Jjq — 30° and s = 1 in this example. Let's 
consider a molecule with a maximum deviation Zi from 
the synchronous molecule. The frequency ujy(t) of this 
molecule can be constructed from its phase (f>{t) as it re- 



volves around the synchronous molecule in longitudinal 
phase-space. Let's now consider an ensemble of molecules 
that is enclosed by this contour, and by a contour that 
is displaced by an infinitesimal value Azi. All these 
molecules experience the same temporal dependence of 
the transverse focusing force, and the transverse trajec- 
tories of the molecules are governed by the equation: 



$ + .;«). ^0. 



(3) 



The transverse phase-space acceptance is easily calcu- 
lated only if ujy{t) is constant. In this case, in which 
LiJy(t) is written as ujy, the longitudinal and transverse 
motions are uncoupled, and in each transverse direction 
the molecules orbit ellipses in transverse phase-space, as 
is shown in FigurefTUl The phase space acceptance {Ay)^^ 
and {Ax)zi in each transverse direction is given by the 
maximum extension ymax — Xmax =1.5 mm from the 
molecular beam axis, and by the maximum transverse 
velocity Vy^max = Vx.max = uiy X ymax that Can be cap- 
tured. The 4D volume {At)zi of the transverse phase- 
space acceptance is then given by 



2n2 



{At)z, = {Ax)zi{Ay)z, = {TTUJyiyniax) ) 



(4) 



If 0Jy{t) is not constant, as is actually the case in a Stark 
decelerator, the molecules experience a transverse fre- 
quency ujy{t) that oscillates between the minimum value 
w™™ and the maximum value uj™^^. The time-averaged 
value of LOy for this molecule is given by: 



(-^). - 



^;{t)dt, 



(5) 



where r is the time it takes the molecule to revolve 
the contour in longitudinal phase-space. The transverse 
phase-space acceptance cannot be calculated anymore, 

, and 



but the values of w,""^, cu'^^'' 



H). 



„ , „,, , ,, ,„„, canneverthe- 

less be used to characterize the transverse phase-space 
acceptance of the ensemble of molecules in three limit- 
ing cases. When the longitudinal oscillation frequency is 
much larger than the transverse oscillation frequency, Uy 

can be taken as ^ (oj'^j . In that case, one obtains the 

4D transverse acceptance for the molecules in this shell 
in longitudinal phase-space. This 4D acceptance can be 
interpreted as the best estimate for the true transverse 
acceptance. When lo™^^ is used, a transverse acceptances 
results that can be interpreted as a strict upper limit for 
the true acceptance. 

The total 6D phase-space acceptance ^(0o) is obtained 
by integrating over all shells with area dAz within the 
separatrix in longitudinal phase-space: 



(6) 



A(0o)= {Ay)zXAx)zAAz 



In Figure [TTJ the longitudinal (2D), transverse (2D) 
and total (6D) phase-space acceptances are shown for 
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the operation modes s = 1 and s — 3. The transverse 
acceptance Ay{(j)Q) is calculated from the total 6D ac- 
ceptance A((f>Q) and the longitudinal acceptance ^^(00) 
via Ay{(l)o) — y^A{(j)o)/A2{(l)o)- For each mode of opera- 
tion three curves are shown; the lower, center, and upper 

curve correspond to the choice of w" 



.mm 



V^R)Z' '''''^ 



for ujy in equation (J4|), respectively. The longitudi- 
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FIG. 11: Prediction for the longitudinal (a), transverse (b) 
and total 6D (c) phase-space acceptance of a Stark decelera- 
tor for the operation modes s = 1 (dashed curves) and s = 3 
(solid curves) as a function of the phase angle (po- For the 
transverse and 6D acceptances, three curves are shown for 
each mode of operation that predict the phase space accep- 
tance in three limiting cases, as explained in the text. The 
center curves describe the best estimate for the acceptance in 
the limit that instabilities in the decelerator can be neglected, 
and have been used in Figure [8] of section IIIII 

nal acceptance ^^((^o) is a factor \/3 smaller for s = 3 
than for s — 1. The transverse acceptance for s ~ 3 
is somewhat larger than for s = 1 and is almost inde- 
pendent of 4>a. Furthermore, the three different curves 
predict a rather similar transverse acceptance for s = 3. 
The three curves for the transverse acceptance for s = 1 
differ much more among each other. The center curves 



that are shown in Figure [TlT c) , and that have been used 
in Figure [S] of section Ifffi predict a rather similar total 
6D phase-space acceptance for s = 1 and s = 3. For al- 
most all phase angles, the lower longitudinal phase-space 
acceptance for s = 3 is compensated for by the larger 
transverse acceptance of s = 3. 



Appendix B 

To obtain a more quantitative value for the thresh- 
old velocity below which excessive focusing occurs in the 
s — 3 mode, we model the motion of a molecule with a 
longitudinal velocity that corresponds to the mean lon- 
gitudinal velocity (vz) of the ensemble. This molecule 
experiences a time-averaged harmonic transverse force 
with a force constant rn(ujz). This force constant is 
deduced from the transverse acceptance Ay{(j)Q) (center 
curve for s = 3 in Figure lllf b) in Appendix A) via 
mloJ^) — mAy{(j)o)/{TTy'^^^)'^, and can be interpreted 
as the weighed average force constant of the ensemble of 
molecules in the packet. 

The transverse force on a molecule in the decelerator 
has periodicity 2sL, and is modeled here by a square 
wave. In a given transverse direction, the molecule re- 
peatedly experiences a constant transverse focusing force 
with harmonic frequency f2j, followed by a section of free 
flight. The sections for focusing and free flight have equal 
length sL, such that Qy = 2(w^). 

For a constant velocity (vz), the resulting equation of 
motion is of the type of Hill's differential equation, and 
the stability of transverse trajectories can be predicted 
with the transfer matrix method. The transformation 
matrices of the coordinate y and velocity Vy for the free 
flight and focusing sections are given by 



Vf 



'Vf 



1 T 
1 



and 



Vf 



'Vf 



coa{nyT) ^ siniflyT) 
-fly sin(51j,r) cos{QyT) 



(7) 



(8) 



respectively, where T — sL/ (vz) and where the sub- 
scripts i and / refer to initial and final coordinates, re- 
spectively. The transformation matrix M for one full 
period is obtained by multiplication of the free flight 
and focusing matrices. The trajectories of the molecules 
through an infinite number of these periods are stable 
if -2 < Tr{M) < 2, where Tr{M) is the trace of the 
transformation matrix [24;]. When Tr{M) is outside 
this range, the amplitude of the transverse motion grows 
without bound. 

In a Stark decelerator, the mean velocity {vz) of the 
packet of molecules is not constant but is gradually 
changed. As long as the velocity change per electric field 
stage is modest compared to the longitudinal velocity, 
the assumptions in the model hold. For phase-angles in 
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the range of 30-40°, ily is found as \/2 • 27r • 720 Hz. For velocities {vz) below 92 m/s. Or, in other words, the 
this value of Qy, the inequality characteristic wavelength A for the transverse oscillation, 

given by A = 27r (vz) /^y, has to be larger than 2.6 • 2sL. 

\2cos{nyT) - nyTsm{nyT)\ < 2 (9) 

implies that transverse trajectories are not stable for 



